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Abstract 

Derivation of the canonical (or Boltzmann) distribution based only on quantum 
dynamics is discussed. Consider a closed system which consists of mutually interacting 
subsystem and heat bath, and assume that the whole system is initially in a pure 
state (which can be far from equilibrium) with small energy fluctuation. Under the 
"hypothesis of equal weights for eigenstates" , we derive the canonical distribution in 
the sense that, at sufficiently large and typical time, the (instantaneous) quantum 
mechanical expectation value of an arbitrary operator of the subsystem is almost equal 
to the desired canonical expectation value. We present a class of examples in which 
the above derivation can be rigorously established without any unproven hypotheses. 

It is often said that the principles of equilibrium statistical physics have not yet been 
justified. It is not clear, however, what statement should be regarded as the ultimate justifi- 
cation. Recalling the astonishingly universal applicability of equilibrium statistical physics, 
it seems likely that there are many independent routes for justification which can be equally 
convincing and important 0]. In the present paper, we concentrate on one of the specific 
scenarios for obtaining canonical distributions from quantum dynamics 0. 

Outline of the work: Let us outline our problem and the main result. We consider an 
isolated quantum mechanical system which consists of a subsystem and a heat bath. The 
subsystem is described by a Hamiltonian which have arbitrary nondegenerate eigenvalues 
£i, . . . ,5„. For convenience we let e^+i > Sj and ei = 0. The heat bath is described by a 
Hamiltonian with the density of states p{B). The inverse temperature of the heat bath 
at energy B is given by the standard formula (3{B) = d\og p{B) / dB . We assume (as usual) 
(3{B) is positive and decreasing in B. The density of states p{B) is arbitrary except for a 
fine structure that we will impose on the spectrum of H^. 

The coupling between the subsystem and the heat bath is given by a special Hamiltonian 
H' which almost conserves the unperturbed energy and whose magnitude is \\H'\\ ~ A. We 
assume ^ A ^ AS, where Ae is the minimum spacing of the energy levels of i^s, and /S.B 
is the maximum spacing of that of H^. These conditions guarantee a weak coupling between 
the subsystem and the bath, as well as macroscopic nature of the bath. The Hamiltonian of 
the whole system is = iJs ® 1b + Is ® H-q + -ff', where Is and 1b are the identity operators 
for the subsystem and the heat bath, respectively. 



^ The first version July 1997. The final version December 1997. Archived as ;ond-mat / 970725^ . To 
appear in Phys. Rev. Lett. 
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Suppose that the whole system is initially in a pure state $(0) which has an energy 
distribution peaked around (but not strictly concentrated at) E. It is possible to treat mixed 
states as well, but such extensions are not essential. For an operator A of the subsystem, we 
denote its quantum mechanical expectation value at time t as 

(A), = (<|.(t),(A®lB)$(t)), (1) 

where (■, ■) stands for the inner product, and $(t) = e~*^*$(0) is the state at time t. Note 
that (■ ■ ■)t is a mixed state on the subsystem. Our main result is the derivation of the 
canonical distribution in the sense that 

for any A (2) 

holds for sufficiently large and typical t, where (^4)^'^° = TTs[Ae~^^^]/TTs[e~^^^] is the 
canonical expectation value. We show that holds for rather general systems under the 
"hypothesis of equal weights for eigenstates." For a special class of models, we prove (Q) 
rigorously without any unproven hypotheses. 

We note the following points about the present derivation of the canonical distribution, 
(i) We do not introduce any probability distributions by hand, (ii) We do not make use of 
the microcanonical distribution, (iii) We do not perform any time averaging, (iv) We do 
not take any limits such as making the bath infinitely large or the coupling infinitesimally 
small, (v) Quantum mechanics seems to play essential roles. 

In the present paper, we describe our main results and basic idea of proofs, leaving 
details to We also briefiy discuss possible extension of the present scenario to more 
general systems. 



Coupling: We diagonalize the (partial) Hamiltonians as Hs'^j = Sj'^j with j = 1, . . . ,n, 
and HB^k = BkTk with k = 1, . . . , N, where "^j, Tk are normalized. We will impose a fine 
structure on the spectrum {Bk} when we discuss our rigorous results. 

When the coupling H' is absent, the total Hamiltonian Hq = ® 1b + Is ® has 
eigenstates Q(j,k) = ® Tfc with eigenvalues f/(j,fc) = £j + B^. We now introduce a new 
index £ = 1, . . . , nN for and U . The index £ is in a one-to-one correspondence with the 
original index (j, k) such that f/^+i > Ui holds for ^ = 1, . . . , nN — 1. We define the coupling 
Hamiltonian H' as 

(G H'Q i) = [ ~ ^'^ ~ f3) 

\ 0, otherwise, 

with a constant A > 0. The Hamiltonian H' describes scattering processes which almost 
conserve the unperturbed energy. 



Eigenstates — main idea: Let $£; be a normalized eigenstate of the total Hamiltonian 
H = Hq + H' with energy E. Expanding it as = J2i=ofi^e, the Schrodinger equation 
E^E = is written as 

E^e = ^{'fe-i + ^e+i) + U^Lpi, (4) 
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with ifo = ifnN+i = 0. This may be regarded as the Schrodinger equation for a single 
quantum mechanical "particle" on a "chain" {1,2,..., nN} under the monotone "potential" 

The eigenstates of (|) with a constant Ue = U are ipi = e^^^^ with k = cos^^[(-E' — U)/X\ 
if \E-U\< A, and ipe = e^"^ (or ifi = (-l)^e^"0 with k = cosh-^[|E - U\/X] ii E - U > X 
{oT E — U < —A). Since our "potential" Ug varies slowly in i, the quasi-classical argument 

suggests that generally (pi takes appreciable values in the "classically accessible region" 
(which consists of i such that \E — Ui\ ^ A), and is negligible outside the region. More 
precisely we expect that, for general E, we can write ^ f{E — Ue) with a {E dependent) 
function f{E) which is nonnegligible only for \E\ ^ A. In terms of the original index, this 
reads 

\^^,,,)\' ^ f[E - {e, + B,)]. (5) 

To see consequences of (^), we take an arbitrary operator A of the subsystem, and denote 
its matrix elements as iA)jj' = {^j,A'^ji). By using (||) and noting that Ae ^ A ^ AB, 
we find 

n N 

j,j'=i k=i 



^j,k \^{j,k)\ 




J2j,k \¥^{j,k)\ 
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^ E,IdBp{B)f{E-e,-B)iA),, 

E,IdBp{B)f{E-6,-B) 
^ p{E - e,){A),, 

^ {A)fCE), (6) 

where the final estimate follows by Taylor expanding p{E — ej) as usual. 

The relation (|^) states that the expectation value in an eigenstate is equal to the desired 
canonical expectation value. This is the key estimate in the present work, and the rest of our 
results follow from relatively general (and standard) arguments. Although we have restricted 
our discussion to the H' of the form (|^), we expect (and can partially prove f^) that the 
property (^ holds for general eigenstates of systems with more general couplings H'. This 
may be called the "hypothesis of equal weights for eigenstates" , from which we may get 
the key estimate (^ and its consequences. 

Eigenstates — rigorous result: We will precisely state the assumption on H^, and de- 
scribe a rigorous estimate corresponding to (|^). We fix an energy unit 6 > (which may be 
much smaller than A), a positive integer R, and positive integers Mi, M2, . . . , M^. We then 
introduce an integer L > Lq (where Lq is a constant), which will be made sufficiently large 
(but finite) to realize the situation where the bath is "large." We require for each r = 1, . . . , R 
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that LMr eigenvalues of are distributed in the interval ((r — 1)5, r5) with an equal spacing 
hr = {LMr)~^S. Thus the density of states of the bath is written as p{r6) = LMj./S. 

With this special Hb, we can partially control the solution of and the sums in (|^) to 
get the following 0. 

Lemma — Consider a system where Hb has the above fine structure, and the coupling 
H' is given by (0). We assume e^+i — ej > 4 A for any j. Let E be an eigenvalue of the 
whole Hamiltonian H such that En + 2X < E < B^i^^ ~ 2A, and be the corresponding 
eigenstates. Then for any operator A of the subsystem, we have 

^EAA®'i-B)^E)-{A)f''\<a\\A\\. (7) 

Here a = 3/3A + 7(£„)2 + cL-^/^^ ^j^^ (3 = (3{E) = d\ogp{E)/dE, 7 = \df3{E - X)/dE\, 
\\A\\ = ma.Xj j' \{A)jj'\, and c is a constant independent of L, E, and A. 

Note that we have a <^ 1 if (i) the coupling is weak (to have /3A <^ 1), (ii) P{E) 
varies slowly (to have 'j{sn)'^ <^ 1), and (iii) the level spacing of the bath is small (to have 
c/,-1/12 ^ Recall that these are the standard assumptions regarded necessary to get the 
canonical distribution. We have established the key estimate (^ rigorously under reasonable 
conditions. 

Long-time average: Let $(0) be the initial state (of the whole system), and expand it as 

m = T.^E^E, (8) 

E 

where the sum runs over the eigenvalues of H. The state at time t is given by ^(t) = 
J2e ^~~^^^1E^E- Then the quantum mechanical expectation value (|l]) can be written as 

{A)t = E e''~^-'''^\^EriE'{<l>E, {A ® 1b)<I>e'). (9) 

E,E' 

Since energy eigenvalues of (^) are nondegenerate, we find that the long-time average of {A)t 
becomes 



E' 

^ T.\iE'\'{Ar,XE'), (10) 

E' 

where F{t) = limr^ooT""^ Jq dtF{t), and we used (||) or (0) to get the final line. 

Further suppose that the initial state $(0) has a small energy fiuctuation in the sense 
that the coefficients je' is nonnegligible only for E' close to some fixed energy E. Then (|10D 
reduces to 

Wt^mTE). (11) 

which states that the long-time average of the quantum mechanical expectation value is 
almost equal to the desired canonical expectation value M. 
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An interesting example of an initial state with small energy fluctuation is 



$(0) = ^„®^afcr,, (12) 

k 



with ak nonnegligible for /c in a finite range. Note that the state (p!2D, restricted onto the 
subsystem, is very far from equilibrium since n is the index for the highest energy state of 
the subsystem. 

Approach to equilibrium: We have seen that, for the initial state with small energy 
fluctuation, the time-independent part in the right-hand side of gives the desired canonical 
expectation value. The time- dependent part of (|^) is a linear combination of many terms 
oscillating in t with different frequencies. It might happen that at sufficiently large and 
typical (fixed) t, these oscillating terms cancel out with each other, and {A)t becomes almost 
identical to its time-independent part (^4)^. 

This naive guess is strengthened by the following simple estimate of the variance. 



E {lE,yiE,{lE.JlE,e'^^'-^'+^'-'^'^HEM\E2){E3\A\E,) 

Ei,E2,E3,E4 

~(EhE\'{E\A\E)] 



E 



< E hE\'hE'\'{E\A\E'){E'\A\E) 

E,E' 

< n^pf max|7£;p, (13) 

E 

where we used the Dirac notation {E\A\E') = {^e,{A ® 1b)^e'), and used the bound 
{E\A'^\E) < n^||y4|p. We also assumed the non- resonance condition for the the relevant 
energy eigenvalues, i.e., whenever Ei — E2 = E4 — E^ ^ holds for Ei such that 7^;. 7^ 0, we 
have El = E4 and E2 = E3. Although we are not able to verify the non- resonance condition 
for a particular given model, it is easily proved that the condition is satisfied for generic 
models [^. When a large number of states equally contribute to the expansion (||), |7£;p and 
hence the right-hand side of ( 1T3|) is very small. This means that {A)t usually takes values 
very close to its average {A)t. 

Following the idea of the Chebyshev's inequality ||10[, this observation can be made 



into a rigorous statement. Instead of writing down the general theorem |^ we present its 
consequence for the special (but interesting) situations with the initial states (|T^). 

Theorem — Suppose that the conditions of the Lemma hold, and the non-resonance 
condition is valid for energy eigenvalues with e„ + 2A < E^ < B^g_^ — 2A. Fix an arbitrary 
energy E with + 3A < < -Bmax "~ 3 A. Take an initial state $(0) of the form (p!2|), and 
assume that au is nonvanishing only for k such that \E — [Bk + en)\ < £^n/2, and satisfies 
(ftfc)^ < c'{enp{E))~^ with an arbitrary constant c' > 1 ||n|. We let a' = 3(3X + 37(£:„)^ + 
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_|_ n^i^Xj e^y/'^ ^ where /3, 7 are the same as in the Lemma. Then there exists a finite 
T > and a subset (i.e., a collection of intervals) Q C [0,T] with the following properties, 
a) For any t G ^, we have 

\{A)t-{A)f^\<a'\\Al (14) 

for any operator A of the subsystem, b) If we denote by yu(^) the total length of the intervals 
in ^, we have 

In simpler words, Q is the collection of "good" time intervals, in which the quantum 
mechanical expectation value of any operator is equal to the canonical expectation value 
within the relative error a' . The new factor A/e„ is small if the coupling is weak. Therefore 
under the physical conditions (i)-(iii) stated below the Lemma, both o' and ?)d{\/ SnY^^ 
in (|15D are small. Then (|l^) says that the "good" intervals essentially cover the whole 
interval [0,T]. We thus have obtained the desired (|^) for most t within the time interval 
< t < T. Recall that the subsystem is far from equilibrium at t = 0. By using only quantum 
dynamics and the assumptions about the initial state, we have established an approach to the 
desired canonical distribution from a highly nonequilibrium state. Although our theorem is 
proved only for systems with particular fine structure and special if', the essential physics 
is contained in the "hypothesis of equal weights for eigenstates" and the non-resonance 
condition (or related weaker conditions). We expect the same scenario to work for much 
more general systems. 

Our result is not strong enough to clarify how {A)t approaches (or deviates from) its 
equilibrium value. Note that we can never expect a perfect decay to equilibrium since {A)t 
is quasiperiodic in t. In a long run, {A)t deviates from the equilibrium value infinitely often. 
Our theorem does guarantee that such deviations are indeed very rare for a weak (but finite) 
coupling. Unfortunately we do not have any meaningful estimate for T. 



Irreversibility: One might question if our theorem implies the existence of an "irre- 
versible" time evolution towards equilibrium. We believe the answer is affirmative, but 
we have to be careful about this delicate issue. We first stress that, exactly as in classical 
cases whether we see irreversibility or not depends on the choice of physical quantities 
that we observe. We should clearly see irreversibility when (there is irreversibility and) we 
observe quantities which have small fiuctuation in the equilibrium (and preferably in the 
initial state too) [|1^ . 



To get an illustrative example, we suppose n is large, and set A to be the projection 
onto the state ^n- A has eigenvalues and 1. In the initial state $(0) of the form (|12|), an 
observation of A does not disturb the state and one gets a definite result 1. We then wait 
for a sufficiently long time (required by the theorem), and once again observe A at time t. 
If this t does not belong to the exceptional "bad" intervals, the theorem guarantees that 
{A)t — {AY^^ = 0{l/n). This does not mean we observe a value of 0{l/n), but means 
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we observe the definite value with probability 1 — 0{l/n) which is essentially 1 for large 
enough n. This is what we mean by (macroscopic) irreversibility |13 . 
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